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The Massive Star Content of NGC 3603 1 

Nicholas W. Melena 2 , Philip Massey 2 , Nidia I. Morrell 3 , Amanda M. Zangari 2 ' 4 

ABSTRACT 

We investigate the massive star content of NGC 3603, the closest known 
giant H II region. We have obtained spectra of 26 stars in the central cluster 
using the Baade 6.5-m telescope (Magellan I). Of these 26 stars, 16 had no 
previous spectroscopy. We also obtained photometry of all of the stars with 
previous or new spectroscopy, primarily using archival HST ACS/HRC images. 
The total number of stars that have been spectroscopically classified in NGC 
3603 now stands at 38. The sample is dominated by very early O type stars 
(03); there are also several (previously identified) H-rich WN+abs stars. We 
derive E(B — V) = 1.39, and find that there is very little variation in reddening 
across the cluster core, in agreement with previous studies. Our spectroscopic 
parallax is consistent with the kinematic distance only if the ratio of total to 
selective extinction is anomalously high within the cluster, as argued by Pandey 
et al. Adopting their reddening, we derive a distance of 7.6 kpc. We discuss the 
various distance estimates to the cluster, and note that although there has been a 
wide range of values in the recent literature (6.3-10.1 kpc) there is actually good 
agreement with the apparent distance modulus of the cluster — the disagreement 
has been the result of the uncertain reddening correction. We construct our H-R 
diagram using the apparent distance modulus with a correction for the slight 
difference in differential reddening from star to star. The resulting H-R diagram 
reveals that the most massive stars are highly coeval, with an age of 1-2 Myr, and 
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of very high masses (120 M ). The three stars with Wolf-Rayet features are the 
most luminous and massive, and are coeval with the non-WRs, in accord with 
what was found in the R136 cluster. There may be a larger age spread (1-4 Myr) 
for the lower mass objects (2O-4OM ). Two supergiants (an OC9.7 I and the Bl 
I star Sher 25) both have an age of about 4 Myr. We compare the stellar content 
of this cluster to that of R136, finding that the number of very high luminosity 
(M hoi < -10) stars is only about l.l-2.4x smaller in NGC 3603. The most 
luminous members in both clusters are H-rich WN+abs stars, basically "Of stars 
on steroids", relatively unevolved stars whose high luminosities results in high 
mass loss rates, and hence spectra that mimic that of evolved WNs. To derive 
an initial mass function for the massive stars in NGC 3603 requires considerably 
more spectroscopy; we estimate from a color-magnitude diagram that less than 
a third of the stars with masses above 2OM have spectral types known. 

Subject headings: HII regions — ISM: individual (NGC 3603) — stars: early-type 
- stars: individual (HD 97950) — stars: Wolf-Rayet 

1. Introduction 

NGC 3603 is our local giant H II region, and our window into the giant H II regions 
seen in other galaxies. Its high luminosity was first recognized by Goss & Radhakrishnan 
(1969). NGC 3603 is located approximately 7-8 kpc from the the Sun in the Carina spiral 
arm. Eisenhauer et al. (1998) estimate that the NGC 3603 cluster has a visual luminosity 
of 6.1 x 1O 5 L , about a factor of five lower than that of the R136 cluster in the heart of 
the giant H II region 30 Doradus in the LMC. The half-light radii are comparable, 4-5 pc 
(Eisenhauer et al. 1998, and references therein). Its relative proximity allows us a unique 
opportunity to study the dense stellar cores that aren't possible in more distant, unresolved 
regions. Its stellar content also serves as a interesting comparison to what is known in R136 
(e.g., Hunter et al. 1996; Massey & Hunter 1998). 

To characterize the massive star population of this high luminosity H II region requires 
both photometry and spectroscopy. The first photometric study of NGC 3603 was by Sher 
(1965), who presented UBV photometry (both photoelectric and photographic) for many of 
the stars in the periphery of the cluster. Sher's (1965) survey was followed by van den Bergh 
(1978)'s photoelectric study, and by Melnick et al. (1989), who used UBV CCD images to 
derive photometry for many of the Sher (1965) stars as well as for some stars located closer 
to to the clusters' core. Melnick et al. (1989) derived a distance of 7.2 kpc, and an age of 
2-3 Myr but with evidence of considerable age spread. Moffat et al. (1994) were able to 
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use HST to obtain F439W (essentially .B-band) photometry of stars even closer to the core, 
although the images were obtained with the original WF/PC camera and suffered from the 
effects of the well-known spherical aberration. Near-IR photometry obtained with adaptive 
optics is discussed by Eisenhauer et al. (1998). Deep JHKL photometry extending down 
to the pre-main-sequence population has been obtained by Stolte et al. (2004), and deep 
optical photometry has been discussed by Sagar et al. (2001) and Sung & Bessell (2004). 
Harayama et al. (2007) further explore the low-mass IMF. Pandey et al. (2000) used multi- 
color photometry to explore the reddening within the cluster, finding that it was anomalously 
high. 

Because of crowding, photometry has been hard enough, but spectroscopy has been all 
but impossible, especially in the central core. HD 97950, the core of NGC 3603, has been 
known to be composite for many years, similar to the situation for R136a in 30 Doradus (see 
Moffat 1983). It has a collective spectral type of around WN6 + 05 (Walborn 1973). Moffat 
(1983) classified 13 objects in the cluster with photographic spectra, finding 11 O-type stars, 
one early B supergiant, plus the unresolved WR core. Some of these spectra were of blends 
that were not resolved. With the Faint Object Spectrograph (FOS) on HST, Drissen et al. 
(1995) were able to classify 14 objects, 11 of which were early O-type stars, including five of 
the earliest class at that time, 03. The spectra of four stars have been modeled: Crowther & 
Dessart (1998) have analyzed spectra of the three Wolf-Rayet stars in the core, and Smartt 
et al. (2002) have studied the B supergiant Sher 25. 

Here matters have stood until the present. Two of the authors (PM and NIM) proposed 
unsuccessfully several times to obtain spectra of stars in the crowded core of NGC 3603 
with the Space Telescope Imaging Spectrograph (STIS) on HST, and after STIS's demise it 
appeared to many of us that the secrets of NGC 3603 would be well-kept for many years 
to come. However, modern ground-based telescopes now are capable of sub-arcsecond im- 
age quality, allowing work from the ground that at one time was only possible from space. 
Here we utilize the good image quality and large aperture of the Baade 6.5-m telescope for 
spectroscopy of individual stars inwards to the dense central core of NGC 3603, classifying 
26 stars, bringing the total number of stars with classification to 38. In addition, we supple- 
ment these spectroscopic data with photometry from recently obtained HST images. These 
data permit a better characterization of the massive star population of this nearby giant 
H II region than has hitherto been possible. In § [2] we will describe the observations and 
reductions. In § [3] we describe our spectral classifications. In § [4] we will use the resulting 
spectral types and photometry to derive new values for the reddening and distances to the 
cluster, and construct an H-R diagram. In § [5] we will summarize and discuss our results. 
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2. Observations and Reductions 
2.1. The Sample 

Previous studies have found that the reddening across NGC 3603 is large (E(B — V) ~ 
1.4) but quite uniform (Moffat 1983), so that the visually brightest stars are likely to have 
the highest visual absolute magnitudes. Since the cluster is dominated by early O-type stars, 
optical photometry cannot tell us much about the physical properties of the stars, as the 
colors are degenerate with effective temperature (see Massey 1998a, 1998b), but does allow 
one to weed out any obvious non-members as was done by Melnick et al. (1989). Thus 
we decided to obtain photometry and spectroscopy for as many of the visually brightest 
members as was practical. To achieve this goal, we identified our sample using HST images 
that had been obtained with the High Resolution Camera (HRC) of the Advanced Camera 
for Surveys (ACS) in connection with a different study of NGC 36ofl We also use these 
images for our photometry (§ 12 . 3j) . 

In Table [1] we list the sample of stars for which we have obtained either spectroscopy 
or photometry. We were able to obtain photometry for all of the stars with spectroscopy 
(both our own spectroscopy and previous), and so Table [T] contains all of the NGC 3603 
stars with optical spectroscopy to the present^. Of these, 26 were obtained by ourselves 
in the present study. We have included a few stars with new photometry for which we 
lack spectroscopy as yet. The coordinates in Table [1] have been determined from the ACS 
images and our own ground-based images, where care has been taken to place these on the 
UCAC2 system (Zacharias et al. 2004). We have retained previous nomenclature for stars 
with previously published spectroscopy or for stars clearly identified in Sher (1965), but 
difficulties in identification from the finding charts published by Melnick et al. (1989) and 
Moffat et al. (1994) have resulted in our presenting new designation for other stars. All stars 
in our sample are identified in Figs. [H and El 



2.2. Spectroscopy 

Optical spectroscopy was obtained by NIM during two nights with the Baade 6.5-m 
Magellan telescope on Las Campanas: 2006 April 12 and 15. The data were taken with the 
Inamori Magellan Areal Camera and Spectrograph (IMACS) in its long camera mode (f/4) 



1 Obtained under program ID 10602, PI=Jesus Mai'z Apellaniz. 

2 This excludes the Moffat (1983) spectra of blobs given letter designations. 
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using the 600 line mm -1 grating. The wavelength coverage was approximately 3600A to 
6700A, with a dispersion of 0.37A pixel -1 , and a spectral resolution of 2.0A. The detector 
for IMACS consists of a mosaic of 8 CCDs. In our case the spectra fell only along one row of 
4 CCDs, and the camera was oriented such that the dispersion crossed all four CCDs, leaving 
three narrow (20A) gaps in our coverage: 4330-4350A (which included H7), 5120-5140A, and 
5920-5940A. A long 0"7 slit was used. 

The spectrograph was usually rotated so that we would obtain two or more stars of 
interest on the slit at the same time. A direct image was usually obtained immediately 
before the spectroscopic exposure, allowing us to be quite certain of the identification of 
stars for which we obtained data, including ones that were coincidentally observed by their 
falling by chance on the slit. Flat field exposures were obtained at the beginning and/or end 
of each night, and HeNeAr comparison arc exposures were taken for each new position. The 
spectroscopic exposures typically consisted of three individual exposures in order to facilitate 
the removal of cosmic rays. Exposure times ranged from 3x250 s to 3x600 s. 

The seeing was good but not spectacular. Using the direct images we obtained adjancent 
to the spectroscopic exposures, we measured an average full width at half maximum of ft'9 
on the first night, and / .'8 on the second. The best images were Qf!67, and the worst were 
I'll. 

Each chip and exposure were reduced separately, and the data were combined at the 
end to produce a single spectrum for each object along the slit. The processing steps were 
the usual ones and were done using IRAF§|. First, the overscan was used to remove the bias 
on each chip, and then the the two-dimensional bias structure was subtracted using the 
average of 10 zero-second exposures. The flat field exposures were divided into the data, 
after normalization. The spectra were extracted using an optimal extraction algorithm, 
after defining the location of stars of interest and selecting the sky background regions 
interactively. The same trace (spatial position as a function of wavelength) and extraction 
apertures were then applied to the comparison arcs, and a wavelength solution was found 
for each aperture. The stellar spectra were then normalized to the continua, combined, and 
the four wavelength regions merged into a single spectrum for each star. 

The typical signal to noise ratio (SNR) is 200-500 per 2A spectral resolution element. 
Such good quality spectra are essential for detecting very weak He I in the earliest O-type 



3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation (NSF). We are grateful to the on-going support of IRAF and the help "desk" maintained by the 
volunteers at |http: / /www.iraf.net| 
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stars. 



2.3. Photometry and Transformation Issues 

Most of our stars were present in the HST ACS/HRC field of view (Fig [2]) and this is 
our primary source for the photometry. The images we analyzed had been taken through 
the F435W and F550M filters, similar to Johnson B and V (but see below). Each image 
consisted of four dithered exposures (offset one from another by a fraction of an arcescond), 
combined into a single "drizzled" image for better sampling of the point-spread function. 
The scale of the final image is 0'.'025 pixel -1 . The total integration time for each image was 
8 s. The images were taken on 2005 December 29. 

For the actual photometry we used the PHOT application of the DAOPHOT package 
of IRAF, adopting an aperture radius of 3.0 pixels (0'.'08). We then determined the aperture 
correction from 3.0 pixels to the ACS "standard" aperture of 0'.'5 using a few isolated stars on 
each frame. We adopted the published aperture correction from 0'.'5 to infinity from Sirianni 
et al. (2005). The counts were corrected for the charge transfer efficiency losses, dependent 
upon the sky background and position on the chip, following the formulation in Pavlovsky 
et al. (2006). We note that the CTE correction is quite significant, amounting to —0.01 to 
—0.07 mag, depending upon the brightness and location of a star. 

Although Sirianni et al. (2005) recommend working in the ACS "native" photometric 
system, we have instead chosen to transform to the standard (Johnson) system. We wish 
to determine the color excesses and distance modulus to the cluster, and to do this we rely 
upon our knowlege of the intrinsic colors and absolute magnitudes of O stars as a function of 
spectral type (see, for example, Conti 1988 and Massey 1998a). These quantities are known 
empirically as the result of observational studies in the standard system rather than through 
models that can simply be recomputed for some alternative bandpass. In other words we 
know that an 07 V star has (B — V)o of —0.32 and My = —4.9 through observations of 
young clusters, not as the result of computations. Thus, in order to compute color excesses 
and derive spectroscopic parallaxes we need to transform something: either we transform 
the intrinsic colors and absolute magnitudes to the ACS system, or we transform the ACS 
photometry to the standard system. We have chosen the latter, since (a) this allows us 
to include stars outside the ACS/HRC field for which standard photometry exists, and (b) 
it allows us to compare our results to others. That said, we need to emphasize that the 
transformations introduce uncertainities well beyond the usual photometric errors. 

Sirianni et al. (2005) list transformations to B for F435W ACS/HRC photometry, but 
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the transformations are based upon unreddened stars, rather than reddened O stars as is 
the case here. Furthermore, they offer no transformation for the F550M filter. So, we 
have determined corrections ourselves as follows. First, we placed the F435W and F550M 
photometry on the so-called "VEGAMAG" system using the revisions to the Sirianni et al. 
(2005) zeropoints listed on the ACS web sitcl- We then corrected the photometry for the 
fact that Vega actually doesn't have zero magnitude, but rather B = V = +0.03 (Bessell et 
al. 1998; see also Mai'z Apellaniz 2006). Using the Kurucz (1992) Atlas 9 models appropriate 
for O-type stars (log g = 5.0 and T eS =40,000, 45,000, and 50,000 K) we then computed 
F550M — V and (F435W — F550M) — (B — V) for various amounts of reddenings, where 
we have adopted the Cardelli et al. (1989) Galactic reddening law. 

For an O star with reddening in the range we expect for NGC 3603 [i.e., E(B — V) — 
1.3 - 1.5], we find F550M — V = -0.11 to -0.12. We adopt F550M — V — -0.11. 
Similarly we find FA35W - B = +0.0 to +0.02, and we adopt F435W - B = +0.01. Thus 
(FA35W - F550M) - (B - V) = +0.120. Jesus Maiz Apellaniz (private communication) 
derives essentially identical corrections for similarly reddened O stars. 

We note that the transformations are sensitive at the 0.02-0.03 mag level to the details 
of the adopted bandpasses. For the ACS bandpasses we synthesized an effective bandpasses 
using SYNPHOT, which includes the wavelength-dependent sensitivities of the entire system 
(telescope + filter + instrument). Historically, knowledge of the standard Johnson band- 
passes has required some "reverse engineering", tested by using the deduced bandpasses 
with model atmospheres to reproduce the observed colors of stars (see Buser & Kurucz 1978 
and Bessell et al. 1998). For B and V SYNPHOT adopts the bandpasses determined by 
Maiz Apellaniz (2006), which are similar to, but not quite identical, to the Bessell (1990) 
versions. Using the Bessell (1990) prescriptions would result in a 0.02 mag shift in both the 
-F435V^ — -B and F550M — V transformations used here. Similarly, if we instead adopted the 
bandpasses determined by Buser & Kurucz (1978), our conversions would differ by 0.03 mag 
for F435W - B and by 0.01 in F550M - V. 

We list the resulting magnitudes and colors in Table [T] as V and B- V. It is worth 
comparing our photometry to that of others, as a check on our transformations. For 9 stars 
in common to Melnick et al. (1989) in our y-band photometry, we exclude two outliers, and 
then find a mean difference of +0.04 mag, in the sense of us minus Melnick et al. (1989). We 
identify only three stars in common for which we have B — V values; for these, we find an 



4 http://www. stsci.edu/hst/acs/analysis/zeropoints. 

5 In contrast, for an unreddened O-type star we find F550M — V = +0.0, FA35W — B = —0.07 and hence 
(F435W - F550N) - (B — V) = -0.07. 
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average difference of —0.02 mag in the same sense. Drissen et al. (1995) list "B" values they 
derived from just averaging the flux in their spectra from 4000A to 4750A. The difference for 



their 13 stars (us minus them) is +0.05 ma, 



We conclude that our photometry is the best 



that we can do unless and until well-calibrated exposures can be made under good seeing 
conditions in filter systems which are a closer approximation to the standard system. 

Some stars fell outside the area covered by the HRC image. Although the same HST 
ACS program included images taken with the Wide Field Camera of ACS, all the stars of 
interest were quite saturated on these archival images. So, instead we used a 20 s V-band 
exposure taken for us by the SMARTS consortium as part of a separate investigation looking 
for eclipsing binaries in NGC 3603. We used three stars in common to set the photometric 
zero-point consistently with the ACS/HRC system. The image was taken on 2006 April 3, 
and the seeing was l'/4, and is used for the finding chart in Fig. [TJ As we lacked a i?-band 
exposure, we adopt the B — V values of Melnick et al. (1989) for those stars, when available. 
The agreement between the SMARTS data and Melnick et al. (1989) is good, with an average 
difference of +0.03 mag, again in the sense of our values minus theirs. The photometry of 
the star Sher 58 clearly differs significantly (AV = 0.58 mag). Their color of this star leads 
them to conclude it is not a member, but their value is inconsistent with the 08 V spectral 
type we find, leading us to suspect that we have identified different stars as Sher 58. The 
star has a close companion, as is seen in Fig. [TJ 



3. Spectral Classifications 

The spectra were classified two ways. First, a qualitative assessment was made by com- 
paring the observed spectra to those illustrated by Walborn & Fitzpatrick (1990). Second, 
a more precise determination was then made quantitatively by measuring the equivalent 
widths (EWs) of He I A4471 and He II A4542 and computing XogW = log(EW He I/EW 
He II). The latter has been calibrated against spectral class by Conti (see Conti & Alschuler 
1971, Conti 1973, Conti & Frost 1977; see summary in Conti 1988), and indeed forms the 
astrophysical basis for the spectral classification of O-type stars. Our experience has shown 
that the visual method works best on spectra with a low or modest SNR, while the latter is 
more accurate for data with a high SNR and adequate spectral resolution. In either scheme 
the primary diagnostic for the spectral subtypes for O stars is the relative strength of He I 
and He II, while the primary luminosity indicators are the strength (emission or absorption) 



6 We exclude the star NGC3603-33, which both they and we find to be marginally resolved on the HST 
images. 
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of the He II A4686 line and the presence and strength of the N III AA4634, 42 emission lines. 
For the earliest O types (02-03.5) we also considered the criteria suggested by Walborn et 
al. (2002), namely the relative strengths of N III and N IV emission, although it has yet to be 
shown whether or not this defines an extension of the temperature sequence (see discussion 
in Massey et al. 2005). We include the new spectral types in Table HJ and illustrate our 
spectra in Fig. [HE 

In a few cases it was clear that our spectral extraction apertures contained some light 
from neighboring objects. If the contamination was judged severe, we did not include the 
star in our study. However, there were a few cases where there was some minor blending, 
and we note such cases. For these, the spectral types may not be as good as for the other 
stars in our sample. 

We were able to classify 26 stars in the cluster, 16 of which were previously unclassified. 
Among these are a number of newly classified 03 and 04 types. We compare our new 
spectral types to those in the literature in Table [21 In general we find very good agreement, 
usually within a single type and luminosity class. Our experience has shown that nebular 
contamination can substantially alter the classification, particularly for the early O stars 
where He I A4387 is not visible, and one must rely solely on He I A4471, which may be 
filled in by nebular emission. In choosing which stars to observe we generally decided not 
to observe the stars that had been previously observed with HST's FOS by Drissen et al. 
(1995). Although the Drissen et al. (1995) sample were observed without the benefits of 
sky subtraction, the FOS observations were obtained with such a small aperture (0'.'25) 
that nebular contamination should not be an issue. However, we did decide to include a 
number of stars that had been previously observed by Moffat (1983), as these were observed 
photographically. Given this, we find the agreement remarkably good, a testament to the 
careful work of the previous studies. That said, there are a few differences we note in 
discussing the stars individually. 

NGC3603-Sh27. Our spectrum of this star (Fig. [5]) is of unusually low signal-to-noise, 
about 130 per 2A spectral resolution element. Nevertheless, it is easy to classify as roughly 
07 V. We measure an EW of 650mA for He I A4471, and 600mA for He II A4542, yielding a 
log W = +0.0 and thus an 07.5 type. The lack of He II A4686 emission makes this a dwarf. 

NGC3603-Sh54. Visually we classified this star as 06 V, with He I A4471 just a bit 
weaker than He II A4542 (Fig. HJ). We measure 310mA and 520mA, respectively, leading to 
a logW = -0.23, also leading to an 06 type (Conti 1988). The strength of He II A4686 
absorption and lack of discernible N III AA4634, 42 emission identifies the star as a dwarf. 

NGC3603-103. The spectrum is shown in Fig. HI This star was slightly blended in 
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our extraction window. We detect extremely weak He I A4471 absorption in our spectrum 
(EW=30mA), typical of the very earliest O type stars (see Massey et al. 2004, 2005). He II 
A4542 has an EW of 660mA, so log W = -1.3. We call this star an 03 V((f)), where the 
luminosity designation denotes a dwarf with weak N III AA4634, 42 emission and He II A4686 
in absorption. 

NGC3603-109. Visually the spectrum of this star appears to be of an 07 V, with He I 
A4471 a bit weaker than He II A4542 (Fig. [5]). We measure EWs of 650mA and 780mA, 
respectively, or \ogW = —0.08. This is consistent with the 07 V designation. The star was 
slightly blended on the slit. 

NGC3603-Sh63. We show the spectrum of this star in Fig. [3j This star has weak He I 
A4471, with a EW of about 180mA, a little large for us to consider the star an 03 (i.e., 
Kudritzki 1980, Simon et al. 1983; see discussion in Massey et al. 2004, 2005). Yet, N IV 
A4058 emission is much stronger than that of N III A4534, 42 emission, which would make 
it an 03 by the criteria enumerated by Walborn et al. (2002). In addition, N V A4603, 19 
absorption is stronger than what we would expect for an 04 star. He II A4542 has an EW 
of 730 mA, and so logW = —0.61, right on the border between 04 and earlier types. In 
addition to N III A4634, 42 emission, He II A 4686 is weak with emission wings, and so we 
call this an 03.5 111(f). Still, N V A4603, 19 appears to be even a bit too strong for this late 
a classification; possibly the star is composite, although here we will treat it as single. The 
star was previously classified by Moffat (1983) as considerably later, and of lower luminosity, 
05.5 V. 

NGC3603-38. The spectrum of this star is shown in Fig. [31 We measure He I A4471 
to have an EW of 120mA, while He II A4542 has an EW of 590mA, leading to an 04 class 
(\ogW = -0.69). There is a little emission at N III AA4634,42 and He II A4686 is weak 
with emission wings, and so we call this an 04 111(f). Previously it was called an 03 V by 
Drissen et al. (1995). The star was slightly blended on our slit, so it is possible that the 
Drissen et al. (1995) type is more accurate. We do see what might be weak He/A4387 and 
A4009, due presumably to the slight blend. 

NGC3603-101. Our visual impression of the spectrum of this star places it in the range 
06-07 V type (Fig. [5]). We measure EWs of 430mA and 580mA for He I A4471 and He II 
A4542, respectively, leading to a log W = —0.13 and an 06.5 type according to Conti (1988). 
The star is a dwarf, judged from the lack of any emission at He II A4686. Very weak emission 
at N III AA 4634, 42 may be present, and we've indicated this by adding a "((f))" description 
to its luminosity class. 

NGC3603-Sh53. Our spectrum of this star is unusually noisy, with a SNR of only 120 
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per 2 A spectral resolution element. Nevertheless, its classification is straight forward as it 
is of mid-0 type, with strong He I A4471 and He II A4542; the former is stronger (Fig. [5]). 
We classify this as an 08.5 V, having measured EWs of 640mA and 380mA, respectively, 
leading to a logW = 0.23. 

NGC3603-104- The spectrum of this star is shown in Fig. [3j He I A4471 is very weak, 
with an EW of 50mA, comparable to that seen in other 03 stars (Kudritzki 1980). The 
EW of He II A4542 is 600mA, leading to a log FT = -1.08. N III AA4634,42 shows weak 
emission, and He II A4686 is weakly in absorption with emission wings, and we call this star 
an 03 111(f). 

NGC3603-Sh56. The spectrum of this star is shown in Fig. [3j This too is a very early 
O-type star, with He I A4471 having an EW of 75mA. Thus we expect this star is of 03 
type. He II A 4542 has an EW of 750mA, leading to logW = —1.00, far more negative 
than the logW = —0.6 used to separate 04s from later type, and thus consistent with our 
assigned type. N III AA4634, 42 is weakly in emission, and He II A4686 has a strong emission 
component (the line appears to be almost P Cygni), suggesting an 03 111(f). However, there 
is a hint of double lines in our spectra, and although we cannot say what the spectral type is 
of the companion, it may be quite early too. We list this spectrum as an 03 III(f)+0? We 
do not include this star when determining the distance modulus in the next section. Moffat 
(1983) classified the star a bit later (04) and of lower luminosity class. 

NGC3603-C. This star shows classic WR features of the WN6 subclass, with strong, 
broad He II A 4686 emission as well as N III AA4634,42 and N IV A4058 emission (Fig. EJ. 
There is absorption superimposed on the Balmer/Pickering lines (H/3 and H5), as well as at 
He II A4542 absorption. We call this WN6+abs, consistent with previous classifications. 

NGC3603-117. The spectrum of this star is shown in Fig. HI The lines in this star are 
quite wide. The EW of He I A4471 is 500mA, while that of He II A4542 is 690 mA, but 
with considerable uncertainty due to the broadness of the lines. With logW = —0.14, the 
Conti (1988) criteria would lean to an 06.5 classification. However, the weakness of He I 
A4387 and general appearance of the spectra suggest a slightly earlier type. We call the star 
an 06 V, with the luminosity class reflecting the lack of emission at N III A4634, 42 and at 
He II A4686. 

NGC3603-Sh25. The spectrum is shown in Fig. |6j This star has long been known to be 
an early-type B supergiant; Moffat (1983) classified it as a B1.5 lab. We would make it just 
slightly earlier, about Bl, based on the relative strengths of Si IV A4089 and Si III A4553. 

NGC3603-Sh64- The spectrum of this star is shown in Fig. HI This is a very early 
O-type star, with He I A4471 just barely discernible on our spectrum; we measure an EW 
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of about 30mA, making the star 03 (or earlier). He II A4542 has an EW of 630mA, making 
logW = —1.3. The luminosity class V, with weak N III AA4634,42 emission but strong 
He II A4686 absorption. We classify it as 03 V((f)). 

NGC3603-102. The spectrum of this star is shown in Fig. [51 The SNR of our spectrum of 
this star is only 100 per 2A spectral resolution element, but fortunately it is easily classified. 
We call it an 08.5 V. He I A4471 has an EW of 950mA, while He II A4542 has an EW of 
530 niA, leading to a logW = 0.25, in agreement with this type. It is a dwarf. 

NGC3603-Sh57. The spectrum of this star is shown in Fig. The EW of He I A4471 
is only 30mA, comparable to what Massey et al. (2004, 2005) found in stars as early as 02. 
He II A4542 has an EW of 540mA, so log W = —1.3. Thus the absorption line spectrum of 
this star suggest it is 03 or earlier. The strength of N III and He II A4686 emission makes 
this a giant. We find that N IV A4058 emission is comparable to N III AA4634,42 emission, 
precluding an 02 classification (Walborn et al. 2002). For giants, Walborn et al. (2002) 
would require this fact to result in an 03.5 111(f) classification, but we are not comfortable 
classifying it this late given the weakness of He I A4471 absorption, and so we call it an 03 
111(f). It was previously called 05 111(f) by Moffat (1983). 

NGC3603-108. The visual impression of this star is that it is roughly 06 V (Fig. H]). 
Measurements of He I A4471 (EW=320mA) and He II A4542 (EW=800 mA) leads to a more 
precise 05.5 (logW = -0.40). There is no emission in N III AA4634,42 or He II A 4686, 
leading to the dwarf luminosity class. The star was slightly blended with a neighbor on the 
slit. We see very weak Mg I A4481 absorption, which is probably due to this blend. 

NGC3603-Shl8. Visually this is a classic 03-04 If star (Fig. EJ), closely resembling 
the spectral standards HD 93129A and HDE 269698 illustrated in Walborn Sz Fitzpatrick 
(1990)0. The strength of N V A4603, 19 is intermediate between the two. N IV A4058 
emission is weaker than N III AA4634,42, which would lean the classification towards an 
03.5 If-04 If type by the criteria listed by Walborn et al. (2002). The EW of He I A4471 
is only 50mA, though, while He II A4542 has an EW of 530 mA, leading to log W = —1.0, 
making this solidly an 03 type using the criteria suggested by Conti (1988). We compromise 
with an 03.5 If type. The strong N III AA4634, 42 and He II A4686 emission leaves no doubt 
to its luminosity class. The star was previously classified as much later, 06 If, by Moffat 
(1983). Nolan Walborn (private communication) reports having independently classified the 
star 04 If from as-yet unpublished data, in agreement with our own type. 

NGC3603-Sh58. Visually this star is roughly 07-08 V (Fig. [5]). We measure nearly equal 



7 Walborn et al. (2002) use HD 93129A now as an example of the 02 If type. 



13 



He I A4471 (EW=640 mA) and He II A4542 (EW=480mA), leading to a log W = +0.13, 
corresponding to an 08 class. He II A4686 is strongly in absorption, and the star is clearly 
a dwarf. The star was somewhat blended on our extraction aperture, and we find that the 
He I lines may be broader than the He II lines, suggesting the resulting type may be a 
composite. We will therefore exclude it when computing the distance modulus in the next 
section. Melnick et al. (1989) listed this star as a non-member (their No. 30) based upon 
their measurement of a very blue color for the star (B — V = 0.47), quite unliked the heavily 
reddened O stars members. However, our photometry gives a color similar to that of the 
other O stars, and a spectral type that confirms membership. 

NGC3603-Sh2l Visually this star is an 06 V (Fig. 0}. We measure an EW of 300mA 
for He I A4471, and an EW of 610 mA for He II A4542. The resulting value log W = -0.31 is 
borderline between an 05.5 and an 06, and we retain the 06 V classification. The strength 
of He II A4686 absorption makes this a dwarf. 

NGC3603-Sh49. Our visual impression of the spectral type of this star is that of an 
08 V (Fig. [5J. We measure an EW of He I A4471 of 570mA, and an EW of He II A4542 of 
560mA, essentially identical, leading us to an 07.5 V type. 

NGC3603-Sh47. This is another early- type O star (Fig. 0}, with weak He I (EW of 
90mA). He II A4542 has an EW of 650 mA, and thus \ogW = -0.86, making the star 04 or 
earlier by Conti (1988). There is only very weak N III A4638,42 emission, and He II A4686 is 
in absorption. The star has previously been called an 04 V by Moffat (1983), and we retain 
this type. 

NGC3603-Sh23. This is a late-type O supergiant (Fig. [6]). The star closely matches the 
spectra of HD 152424 and HD 104565 shown by Walborn Sz Fitzpatrick (1990), and we thus 
classify the star as an OC9.7 la, where the "C" denotes the excessively strong CIII A4650 
line. This is in substantial agreement with the 09.5 lab type found by Moffat (1983). 

NGC3603-Sh22. This is clearly a very early O-type star (Fig. [3]). Despite a signal-to- 
noise ratio of >350 per spectral resolution element, we detect no He I A4471. The EW must 
be <20mA. N IV A4058 and N III AA4634,42 emission are both weak, but unfortunately 
the N IV line also coincides with some bad pixels. We call the star an 03 111(f), although 
it could be an 02 by the Walborn et al. (2002) criteria. The giant luminosity class comes 
about from the modest emission at He II A4686. The star had been called an 05 V (f) by 
Moffat (1983). 

NGC3603-Sh 21 . The spectrum of this star is shown in Fig. [51 Our initial impression of 
the spectrum of this star is that it is roughly of type 06 V, with He I a bit weaker than He II. 
We measure EWs of 390 mA and 700 mA, respectively, leading to a value logW 7 "' = —0.25, 
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consistent with the 06 class. 

NGC3603-Sh 19. The spectrum of this star is shown in Fig. |H This is another very 
early O star. Our spectrum has a fairly low SNR (200 per 2A resolution element) and we 
detect no He I. We call this an 03 V((f)). 

4. Results 

We are now prepared to derive values for the reddening, distance, and age of the 
NGC 3603 spectroscopic sample. We begin with the reddenings. In Table [3] we include 
the values of E(B — V) we derive, using the intrinsic colors of Table 3 from Massey (1998a), 
and the observed colors from Table [H 

We derive an average reddening E(B — V) = 1.394 ± 0.012 (standard deviation of the 
mean, hereafter "s.d.m."). The median is 1.39. This is very similar to Moffat (1983)'s 
finding of an average value of E(B — V) — 1.44. Sung & Bessell (2004) find a somewhat 
lower value, E(B — V) — 1.25, in the core of the cluster, but argue that it increases to larger 
values at greater distances. Moffat (1983) noted that there was very little spread in the 
reddening among the cluster stars — his sample showed a dispersion of only 0.09 mag. We 
find an even smaller dispersion, 0.06 mag. We are thus reassured that the reddening is well 
determined to the O stars in NGC 3603, and agree with the conclusion of Moffat (1983) that 
there is very little variation in reddening across this cluster. We adopt an average value of 
E(B -V) = 1.39. 

It is interesting to note that the Bl supergiant, Sher 25, has a color excess that is 
considerably larger than average. This supergiant has been compared to the progenitor of 
SN 1987A (see Smith 2007 and Smartt et al. 2002), and is known to show circumstellar 
material (Brandner et al. 1997). 

Next, let us derive a distance modulus to the cluster by adopting an absolute magnitude 
for each star based upon the spectral type and luminosity class. We adopt the values of Conti 
et al. (1983) for the O stars, interpolating as needed; for the B supergiant (Sher 25) we adopt 
My = —6.5 from Humphreys Sz McElroy (1984). If we make no correction for the reddening 
to the cluster, we derive an apparent distance modulus of 19.12±0.09 (s.d.m.). The standard 
deviation of this sample is 0.5 mag, which is what we expect, given that that is also the typical 
scatter in My for a given spectral type and luminosity class (Conti 1988). The median is 
19.00. If we exclude the giants and supergiants (as their absolute magnitudes might cover 
a larger range), we determine similar values: the average is 19.21 ± 0.12 (s.d.m.), with a 
standard deviation of the sample of 0.6 mag. The median is 19.07. We therefore adopt an 
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apparent distance modulus to the cluster of 19.1 ± 0.1. 

If the extinction were normal (Ry = 3.1) this would then correspond to a true distance 
modulus of 14.8, or 9.1 kpc. However, Pandey et al. (2000) have investigated the reddening 
towards this cluster, and conclude that the reddening is anomalously high within, with a ratio 
of total to selective extinction Ry = 4.3, a value that is more typical of dense environments 
(see, for example, Whittet 2003). They correct for reddening assuming an E(B — V) = 1.1 
for foreground (with R v = 3.1), and Ry = 4.3 for the color excess above this value; i.e., 
Ay = 3.41 — 4.3 [E(B — V) — 1.1]. If we were to make this correction, then the true distance 
modulus we would obtain would be 14.4, or 7.6 kpc. 

The physical distance to NGC 3603 is poorly determined. There are three methods that 
have been employed: main-sequence fitting, spectroscopic parallaxes, and kinematic distance 
determinations based upon rotation models of the Milky Way. We summarize the results 
in Tables H] and [5j Main-sequence fitting (Melnick et al. 1989, Pandey et al. 2000), Sung & 
Bessell 2004) is probably the least certain of these for such a young cluster, given that the 
theoretical zero-age main-sequence is nearly vertical in the color-magnitude diagram. This 
is particularly true in the V — I plane used by Sung & Bessell (2004); see their Fig. 9. This 
problem was recognized by van den Bergh (1978), who obtained UBV photometry of cluster 
members, but noted that such data was not sufficient to determine a reliable distance. (He 
did, however, determine a kinematic distance, as discussed below.) See also the discussion 
in Sagar et al. (2001). The spectroscopic parallaxes (Moffat 1983, Crowther & Dessart 1998, 
Sagar et al. 2001) should be more reliable. Our sample of stars is considerably larger than 
those previously employed. Of course, a key issue in deriving the physical distance is how 
reddening is treated. In Table H] we include the apparent distance moduli as well as the 
derived physical distances. We can see that most modern studies derive a similar value for 
the apparent distance modulus (i.e., 18.6-19.1), despite the extremely large range in physical 
distances derived (6.3 kpc to 10.1 kpc). 

We can compare our distance to the kinematic distances (Table |5]). Goss & Radhakrish- 
nan (1969) find a distance of 8.4 kpc, but that was based upon the old assumption that the 
distance from the Galactic center to the sun is 10 kpc. Correcting to the more modern value 
Rq = 8.5 kpc would result in a 7.1 kpc distance. Van den Bergh (1978) also derives a similar 
kinematic distance, 6.8 kpc, when corrected to Rq = 8.5 kpc. De Pree et al. (1999) quote a 
value of 6.1 ±0.6 kpc, but one of the co-authors, Miller Goss (private communication), writes 
that this value should be revised to 7.0 kpc. Russeil (2003) finds a slightly larger kinematic 
distance, 7.9 kpc. While these span a range, we can conclude that our apparent distance 
modulus can be reconciled with the kinematic distance moduli in the literature only if the 



-16- 



extinction is anomalous within the cluster, as Pandey et al. (2000) arguej. 

For constructing the H-R diagram we have chosen to adopt our apparent distance 
modulus, and then correct individual stars for differences in the average reddening by 
4.3 x AE(B — V), where AE(B — V) is the difference between the average color excess 
1.39 and that for an individual star. Since the variation in the reddening is small, using 3.1 
instead would lead to at most a difference of 0.25 mag in the final value. We include the 
absolute visual magnitudes obtained by this method in Table [3J In most cases the agreement 
between the My (based on the spectral type and luminosity class) and that computed from 
the average apparent distance modulus (corrected for reddening above or below the cluster 
average) is good, within the 0.5 mag scatter we expect (Conti 1988). 

Finally, let us use the recent effective temperature scale of Massey et al. (2005) to 
place these stars on the H-R diagram. The new scale includes the effects of line- and wind- 
blanketing, which significantly lowers the effective temperatures for Galactic stars. For the 
Bl I star (Sh 25) we adopted an effective temperature of 22000 K, consistent with the 
modeling of Smartt et al. (2002) and the recent findings of Crowther et al. (2006). For 
that star we have adopted a bolometric correction of -2.0 mag, also consistent with these 
studies. We have included the three WR stars based upon the modeling of Crowther & 
Dessart (1998), where we have increased the luminosity by 0.12 dex for consistency with 
the larger apparent distance modulus found here. We show the resulting H-R diagram in 
Fig. [7J We compare these to the Schaller et al. (1992) evolutionary tracks, which correspond 
to "Galactic" (solar) metallicity (z = 0.020), which is appropriate for NGC 3603 (Smartt 
et al. 2001; Peimbert et al. 2007; Lebouteiller et al. 2007). The black solid lines denote the 
isochrones at 1 Myr intervals from an age of 1 Myr to 6 Myfl 

First, we see that the masses range above 12OM . We see that the most massive stars 
are the stars with Wolf-Rayet features analyzed by Crowther & Dessart (1998). In this 



8 Smartt et al. (2002) reject Ry = 4.3 for Sher 25 based on the argument that applying it leads to too high 
a luminosity for this star, and instead adopt Ry = 3.7 on somewhat arbitrary grounds. However, they had 
(mis)applied Ry — 4.3 to the entire line of sight reddening (E(B — V) = 1.60), rather than applying Ry = 3.1 
to the foreground reddening (E(B - V) = 1.1) and Ry = 4.3 to the remainder (AE(B - V) = 1.60 - 1.10). 
Thus, using the Pandey et al. (2000) reddening leads to Ay — 5.56 mag, not 6.88 mag, and actually smaller 
than the value Smartt et al. (2002) adopt {Ay = 3.7x1.60=5.92 mag). 

9 We have used the older Geneva evolutionary tracks of Schaller et al. (1992), rather than the newer ones 
of Meynet & Maeder (2003), simply for convenience, as Georges Meynet had kindly made available software 
to compute isochrones from the older tracks. The primary difference with the newer tracks is the inclusion 
of rotation. While rotation significantly alters the tracks of low-metallicity stars, such as those in the SMC 
and the LMC, there is much less of an effect at Galactic metallicities; see Meynet & Maeder (2000). We 
have also truncated the tracks at the start of the WR stage just for clarity. 
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context it is worth remembering that in the R136 cluster there were several H-rich WN6 
stars which had absorption lines; Massey & Hunter (1998) concluded that these were not 
evolved objects, but rather "Of stars on steroids", i.e., stars whose masses were so high, and 
which were so luminous, that their winds were so strong that their spectra simply resembled 
WR stars in having strong emission lines. The modeling of Crowther & Dessart (1998) bear 
this out: the WRs are coeval with the rest of the cluster, and are simply slightly more 
luminous and massive. (Note that for simplicity we have truncated the evolutionary tracks 
prior to the Wolf-Rayet stage.) The most massive non-WR stars in R136 are more massive 
than we see here, but as discussed below, that could be a metallicity effect. In any event, 
the unevolved stellar population of NGC 3603 is certainly quite massive. 

Secondly, the ages of the highest mass stars are roughly 1-2 Myr, with very little dis- 
persion. Sung & Bessell (2004) obtained a similar result. Stars of somewhat lower mass 
(< 40M Q ) show a larger dispersion in age, with ages ranging from 1-4 Myr. The evolved 
B-type supergiant Sh 25 and the OC9.7 star Sh 23 have an age of about 4 Myr. This is 
similar to the age found for stars in the outer regions of the cluster by Sung & Bessell (2004). 

It is premature to attempt to derive an initial mass function for the high mass stars 
of NGC 3603, as photometry alone cannot provide a sufficiently accurate discriminant of 
bolometric luminosity (and hence mass) at such high effective temperatures. (See discussion 
in Massey 1998a, 1998b.) Sagar et al. 2001 and Sung & Bessell 2004 do use their deep 
photometry to compute IMFs for the intermediate-mass stars, but for these photometry 
does provide a good discriminant. In Fig. [8] we show the color magnitude diagram of the 
cluster, where we have indicated the stars for which there are spectral types by open circles. 
We have cut the diagram at V = 15.6, corresponding to a ZAMS 20M© star. (Such a 
star would have My = 3.5 and would correspond to a late O-type near the ZAMS.) The 
photometry for the general field (shown in black) comes from Melnick et al. (1989), while 
that of the central field (shown in red) comes from the HRC image. We see that although 
we have made significant progress in the spectroscopy of stars in NGC 3603, much work 
remains to be done. In the area outside the core, there are spectra for about 20% of the 
interesting stars(l/ < 15.6 with colors indicating likely membership), while in the central 
portion we have spectra for about a third of the stars with V < 15.6. (Note that all of the 
stars in the core have colors indicative of membership). In all, spectra of another 90 stars 
would be needed for a complete census down to 2OM . In the central region, observing the 
majority of the remaining stars will be difficult due to crowding, but not impossible. We 
hope to undertake such work during the next observing season. 
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5. Discussion and Summary 

It is worth comparing what we know now about the massive star content of NGC 3603 
with that known about the R136 cluster at the heart of the giant H II region 30 Doradus in 
the LMC. Spectroscopy of the stars in R136 is similarly complete (Massey & Hunter 1998). 
R136 contains an even greater wealth of 03 stars, and stars extending up to 150M Q , about 
the point where the IMF peters out to a single star. In the H-R diagram shown by Massey & 
Hunter (1998), most of the massive stars appear to be strictly coeval, with an age between 1 
and 2 Myr, just as we find here for NGC 3603. For stars with masses below about 4OM the 
placement of stars in the H-R diagram of R136 was constrained only by photometry, and so 
there is an (apparent) spread in ages from 0-6 Myr, although the actual degree of coevality 
may be higher. There is a single B-type supergiant in their H-R diagram with a mass of 
about 2OM , and an age of 10 Myr. Maybe that star is an interloper from the field of the 
LMC, or it could be that in both NGC 3603 and in R136 we simply see that a couple of high 
mass stars formed a bit earlier than the majority of stars in the cluster. 

We note that both of these clusters contain H-rich WN+abs stars. Massey & Hunter 
(1998) argued that these were unevolved (i.e., core H-burning) stars of high luminosity and 
mass, whose spectra mimicked that of WN stars given the strong stellar winds expected from 
high luminosity. Given the ages of the high mass stars in NGC 3603 or in R136, it is not 
possible for the WR stars to be evolved (core He-burning) objects, not if they are coeval with 
the rest of the massive stars. The physical properties of these stars determined by modeling 
by Crowther & Dessart (1998) is consistent with thig^l At Galactic metallicity we would 
expect the onset of such features to happen at lower mass than in the LMC, so these "Of 
stars on steroids" in NGC 3603 may be less massive than in the R136 cluster. 

We can compare the stellar content straightforwardly. Let us simply count the number 
of stars with known spectroscopy brighter than a certain bolometric luminosity, Mboi ~ —10. 
(For these, our spectroscopy is mostly complete; (V ~13-13.5; see Fig. [8]). From Table [3] we 
find 9 stars listed in NGC 3603 — 12, if we then include the WN+abs stars in the tally. There 
are perhaps another 6 stars that have not been observed spectroscopically that could be as 
bolometrically luminous, so the total is 12-18. In R136 we find 20-29 such stars, depending 
upon the adopted temperature scale (i.e, Table 3 of Massey & Hunter 1998). So, while R136 
is richer in massive stars than NGC 3603, it is only by a factor of 1.1 to 2.4. Moffat et al. 
(1994) have also argued for the similarity of the central clusters, suggesting that the primary 
difference is that NGC 3603 lacks a surrounding massive halo of cluster stars. 



10 These stars were found to be chemically evolved, in the sense of showing enhanced CNO products at the 
surface, but this is still consistent with these being H-burning objects. 
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In summary, we have obtained spectra of 26 stars in the NGC 3603 cluster, 16 of which 
have no previous spectroscopy. That brings the total number of stars with spectral types to 
38 (Tabled]). In addition, we provide identification and ACS/HRC photometry of another 
12 stars in the central core for which spectroscopy would be desirable. Our spectroscopic 
sample includes many stars of type 03. We find an average reddening E(B — V) = 1.39, with 
very little scatter (sample standard deviation of 0.05 mag), indicative of very little variation 
in reddening in the core. Our spectroscopic parallax for the cluster can be reconciled with 
the kinematic distance of the cluster if we adopt the reddening proposed by Pandey et 
al. (2000); we then derive a distance of 7.6 kpc. We emphasize that although there has 
been a very large range of physical distances derived for this cluster in the past 10 years 
(6.3-10.1 kpc), there is excellent agreement in the apparent distance modulus of the cluster 
(18.6-19.1 mag). The disagreement in the distances are really based upon exactly how to 
correct for reddening. We side step the issue by using the apparent distance modulus (plus 
a modest correction for differential reddening) to construct the H-R diagram. It reveals a 
mostly coeval population of massive stars extending beyond 12OM , with ages of 1-2 Myr. 
The most massive and luminous stars are the H-rich WN stars, as expected by analogy with 
R136 (Massey & Hunter 1998). Some stars of 2O-4OM may show a larger age spread, up to 
4 Myr, with the OC9.7 and B-type supergiants having an age of 4 Myr. 

Additional spectroscopy of stars is underway by A. F. J. Moffat (private communica- 
tion), which will further help refine our view of this interesting cluster. While further HST 
STIS spectroscopy of the most crowded members would be highly desirable, it is worth not- 
ing how much can now be accomplished from the ground using the best modern telescopes. 
More such work is needed for deriving an initial mass function for the massive stars in this 
cluster, as less than a third of the stars with masses above 20M Q have spectroscopy. We 
hope to help improve this situation in the coming observing season ourselves. 

We are grateful to David Osip for allowing us to observe during some Magellan engineer- 
ing time shortly in advance of our scheduled night. Jesus Maiz Apellaniz was kind enough 
to allow us access to one of the ACS/HRC frames prior to it becoming public in order to 
prepare for the observing run; we also gratefully acknowledge his calling our attention to an 
error we made in the photometric transformations in an earlier version of the manuscript. 
Jennifer Mack provided guidance in our use of SYNPHOT. We also benefited from corre- 
spondence with Geoff Clayton, Miller Goss, Tony Moffat, Brian Skiff, and Nolan Walborn, 
as well as from comments on an early draft by Deidre Hunter. Suggestions by the referee, 
Stephen Smartt, were useful in improving the paper. NWM and PM work was supported 
through the National Science Foundation (NSF) AST-0604569. AMZ's work was supported 
through an NSF REU grant, AST-0453611. 
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Fig. 1. — . Finding chart for stars in our sample from the periphery of NGC 3603. The field 
of view is about l'O on a side. North is up, and east is to the left. The chart is made from a 
V-band image taken with the CTIO Yale 1.0-m telescope. The circles are 4'.'6 in diameter, or 
0.17 pc at an assumed distance of 7.6 kpc. The square denotes the approximate area shown 
in Fig. El 
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Fig. 2. — . Finding chart for stars in our sample from the center of NGC 3603. The field 
of view is 20"on a side, and roughly corresponds to the square outline shown in Fig. CD The 
image has been rotated so that north is up and east is to the left. The chart was made 
from a V-band ACS/HRC image. The circles are 0'.'4 in diameter, or 0.02 pc at an assumed 
distance of 7.6 kpc. 
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Fig. 3. — Normalized spectra of NGC 3603 supergiants and giants. Only the data in the 
blue (MK classification) region are shown. The interstellar (is) H and K Ca II lines and 
diffuse interstellar bands (dib) are marked, along with the prominent stellar features. The 
displacement between the various spectra is 0.4 times the continuum level. 
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Fig. 4. — Same as Figure [3] for 03-06 dwarfs. 
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Fig. 5. — Same as Figure [3] for 06-08.5 dwarfs, except that the displacement between the 
spectra is 0.6 times the continuum level. 



Fig. 6. — Same as Figure [3] for the two "late-type" supergiants in our sample, and OC9.7 
lab star (Sh 23) and a Bl lab star (Sh 25). The bluest part of the spectrum in the latter 
star has been heavily smoothed. The displacement between the two spectra is 0.6 times the 
continuum level. 
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Fig. 7. — The H-R diagram for NGC 3603. Solid points show the location of the stars 
listed in Table [3j The non-rotating evolutionary z = 0.020 tracks of Schaller et al. (1992) 
are shown, with the initial masses labeled, in various colors. Isochrones are shown as black 
curves for 1, 2, 3, 4, 5, and 6 Myr. We have labeled the locations of OC9.7 la star Sh23, and 
the Bl lab star Sh 25. The three asterisks denote the location of the three H-rich Wolf-Rayet 
stars, with the physical properties taken from Crowther & Dessart (1998). 
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Fig. 8. — The color-magnitude diagram for NGC 3603. Open circles denote the stars for 
which there are spectral types; crosses denote the stars for which there are only photometry. 
Red indicates the stars in the central region (i.e., that covered in Fig. EJ), and are based upon 
point spread function fitting to stars on the HRC images. (This photometry will be fully 
presented in a later paper.) Black indicates photometry of stars outside the central region, 
from Melnick et al. (1989). 



Table 1. NGC 3603 Stars with New Data 



Star a 02000 ^2000 V B - V Spectral Type Rcf 



Sh27 


11 


15 


03 


921 


-61 


15 


23 


06 


15 


04 


1 


05 b 


07.5 V 


1 


Sh54 


11 


15 


03 


976 


-61 


15 


35 


72 


14 


57 


1 


08 b 


06 V 


1 


103 c 


11 


15 


06 


237 


-61 


15 


36 


58 


13 


09 





99 


03 V((f)) 


1 


128 


11 


15 


06 


422 


-61 


15 


36 


42 


14 


76 





98 






109 c 


11 


15 


06 


588 


-61 


15 


40 


40 


13 


85 





99 


07 V 


1 


Sh50 


11 


15 


06 


615 


-61 


15 


23 


81 


14 


70 


1 


14 






129 


11 


15 


06 


685 


-61 


15 


38 


85 


14 


75 





99 
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11 


15 


06 


726 


-61 


15 


39 


26 


14 


96 





99 






111 


11 


15 


06 


741 


-61 


15 


35 


64 


13 


87 


1 


00 






124 


11 


15 


06 


815 


-61 


15 


33 


95 


14 


54 


1 


02 






Sh63 


11 


15 


06 


847 


-61 


15 


44 


69 


13 


41 


1 


04 


03.5 111(f) 


1 


38 c 


11 


15 


06 


917 


-61 


15 


36 


60 


13 


21 





97 


04 111(f) 


1 


45 


11 


15 


06 


944 


-61 


15 


38 


30 


14 


14 


1 


00 


08 V-III 


2 


37 


11 


15 


07 


003 


-61 


15 


37 


45 


14 


16 





97 


06.5+OB? 


2 


41 


11 


15 


07 


046 


-61 


15 


38 


95 


14 


24 


1 


00 


04 V 


2 


42 


11 


15 


07 


060 


-61 


15 


39 


30 


12 


99 


1 


03 


03 III(f*) 


2 


101 


11 


15 


07 


068 


-61 


15 


45 


32 


14 


02 


1 


08 


06.5 V((f)) 


1 


36 


11 


15 


07 


078 


-61 


15 


37 


56 


14 


52 





99 


06 V 


2 


40 


11 


15 


07 


124 


-61 


15 


39 


09 


13 


33 


1 


05 


03 V 


2 


Sh53 


11 


15 


07 


148 


-61 


15 


54 


86 


14 


47 


1 


12 b 


08.5 V 


1 


120 


11 


15 


07 


211 


-61 


15 


41 


50 


14 


35 


1 


03 






104 


11 


15 


07 


279 


-61 


15 


34 


94 


13 


02 


1 


06 


03 111(f) 


1 


Al 


11 


15 


07 


305 


-61 


15 


38 


43 


11 


18 


1 


03 


WN6+abs 


2 


A2 


11 


15 


07 


313 


-61 


15 


38 


79 


12 


53 


1 


04 


03 V 


2 


A3 


11 


15 


07 


352 


-61 


15 


38 


46 


13 


09 


1 


04 


03 III(f*) 


2 


33 d 


11 


15 


07 


363 


-61 


15 


39 


54 


13 


69 


1 


03 


05 V+OB? 


2 


B 


11 


15 


07 


411 


-61 


15 


38 


58 


11 


33 


1 


01 


WN6+abs 


2 


Sh56 


11 


15 


07 


498 


-61 


15 


46 


35 


13 


48 


1 


08 


03 III(f)+0? 


1 


C 


11 


15 


07 


589 


-61 


15 


38 


00 


11 


89 


1 


05 


WN6+abs 


1 


139 


11 


15 


07 


591 


-61 


15 


32 


64 


14 


88 


1 


08 
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11 


15 


07 


593 


-61 


15 


35 


96 


14 


10 


1 


04 






117 


11 


15 


07 


623 


-61 


15 


30 


24 


14 


17 


1 


10 


06 V 


1 


125 


11 


15 


07 


630 


-61 


15 


31 


35 


14 


50 


1 


06 






Sh25 


11 


15 


07 


649 


-61 


15 


17 


59 


12 


23 


1 


42 b 


Bl lab 


1 


141 


11 


15 


07 


692 


-61 


15 


34 


65 


15 


02 


1 


04 






Sh64 


11 


15 


07 


822 


-61 


15 


27 


93 


13 


58 


1 


15 


03 V((f)) 


1 


16 


11 


15 


07 


825 


-61 


15 


37 


84 


13 


53 


1 


08 


03 V 


2 
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Table 1 — Continued 



Star a 




«2000 






C>2000 




V 


B 


- V 


Spectral Type 


Rcf 


102 


11 


15 


08 


.073 


-61 


15 


34. 


.24 


15.32 


1 


.10 


08.5 V 


1 


135 


11 


15 


08 


.106 


-61 


15 


40. 


.30 


14.75 


1 


.06 






Sh57 


11 


15 


08 


.198 


-61 


15 


47. 


.33 


13.28 


1 


,06 b 


03 111(f) 


1 


108 c 


11 


15 


08 


.513 


-61 


15 


38. 


.30 


13.71 


1 


.10 


05.5 V 


1 


Shl8 


11 


15 


08 


.712 


-61 


15 


59. 


.95 


12.65 


1 


,21 b 


03.5 If 


1 


Sh58 c 


11 


15 


08 


.699 


-61 


15 


44. 


.49 


14.24 






08 V° 


1 


Sh24 


11 


15 


08 


.905 


-61 


15 


27. 


.32 


14.27 


1 


,ll b 


06 V 


1 


Sh49 


11 


15 


09 


.129 


-61 


15 


33 


,18 


14.67 


1 


,10 b 


07.5 V 


1 


Sh47 


11 


15 


09 


.353 


-61 


16 


02. 


,07 


12.72 


1 


,14 b 


04 V 


1 


Sh23 


11 


15 


09 


.849 


-61 


15 


30 


,48 


12.70 


1 


,14 b 


OC9.7 la 


1 


Sh22 


11 


15 


10 


.071 


-61 


15 


38. 


,01 


13.21 


1 


,06 b 


03 111(f) 


1 


Sh21 


11 


15 


11 


.074 


-61 


15 


36. 


,85 


14.75 


1 


,12 b 


06 V((f)) 


1 


Shl9 


11 


15 


11 


.317 


-61 


15 


55. 


,63 


13.65 


1 


,02 b 


03 V((f)) 


1 



Note. — Units of right ascension are hours, minutes, and seconds, and units 
of declination are degrees, arcminutes, and arcseconds. 

Nomenclature: stars with "Sh" designations are from Sher 1965. Stars with 
two digits are from Moffat et al. 1994, and are re-identified in Drissen et al. 
1995 as "HST-nn" . The "A" , "B" , and "C" components in the central core are 
from Drissen et al. 1995. The three-digit numbers are from the present paper. 

h B - V values from Melnick et al. 1989. 

c Slightly blended on slit. 

d Blended on ACS image. 

e Composite? 

References. — For spectral types: (1) This paper; (2) Drissen et al. 1995. 
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Table 2. NGC 3603 Stars with Newly Found Spectral Types 



Star 


New 


Literature 






Spectral Type 


Spectral Type 


Ref 


Sh27 


07.5 V 






Sh54 


06 V 






103 a 


03 V((f)) 






109 a 


07 V 






Sh63 


03.5 111(f) 


05.5 V 


1 


38 a 


04 111(f) 


03 V 


2 


101 


06.5 V((f)) 






Sh53 


08.5 V 






104 


03 111(f) 






Sh56 


03 III(f)+0 ? 


04 V(f) 


1 


C 


WN6 + abs 


WN6 + abs 


2 


117 


06 V 






Sh25 


Bl lab 


B1.5 lab 


1 


Sh64 


03 V((f)) 






102 


08.5 V 






Sh57 


03 111(f) 


05 111(f) 


1 


108 a 


05.5 V 






Shl8 


03.5 If 


06 If 


1 


Sh58 


08 V a 






Sh24 


06 V ••• 






Sh49 


07.5 V 






Sh47 


04 V 


04 V 


1 


Sh23 


OC9.7 la 


09.5 lab 


1 


Sh22 


03 111(f) 


05 V(f) 


1 


Sh21 


06 V((f)) 






Shl9 


03 V((f)) 







References. - References for spectral types 
from the literature: (1) Moffat 1983; (2) Drissen 
et al. 1995. 

a Composite? 
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Table 3. Reddening, Distances, Absolute Magnitudes 



Star 


V 


B-V 


Spectral 


E(B - V) 


Adopted 


V - My 


Computed 


l0gT 0ff 


M bol 








Type 








M v b 






Sn27 


15.04 


1 A ET 

1.05 


(J7.5 V 


1 07 

1.37 


CT A 

-5.0 


OA A A 

20.04 


A A 

-4.0 


A E A 1 

4.541 


T O 

-7.3 


Sii54 


14.57 


1 AO 

1.08 


(J5 V 


1 A A 

1.41) 


-5.1 


1 A C7 

19.67 


-4.6 


A E D O 

4.583 


O O 

-8.2 


1 AO 

103 


1 nn 

13.09 


A A A 

0.99 


03 V((fj) 


1.31 


-5.4 


IO A A 

18.49 


-5.7 


4.667 


A A 

-9.9 


109 


13.85 


A A A 

0.99 


07 V 


1.31 


A A 

-4.9 


18.75 


A A 

-4.9 


4.556 


O A 

-8.4 


Sn63 


13.41 


1.04 


r\ r ttt/'-C\ 

03. 5 111(1) 


1.37 


-5.9 


19.31 


-5.6 


4.653 


-9.7 


OO 

38 


13.21 


A AT 

0.97 


TTTZ-f^ 
(J4 lll(r) 


1 OA 

1.30 


a a 
-5.0 


1 A Ol 

19. 21 


-5.5 


A £Z A O 

4.543 


a e 

-9.5 


A C 

45 


1 A 1 A 

14.14 


1 AA 
1.00 


C\Q \T TTT 

(Jo V-lll 


1 00 
1.62 


-4.0 


1 n a 

10.94 


A T 

-4.7 


/I coo 


T A 

-7.9 


A 1 

41 


1 A O A 

14.24 


1.00 


04 V 


1 00 
1.62 


-5.8 


OA A A 

20.04 


-4.0 


/I £Z A O 

4.543 


a 

-8.6 


42 


12.99 


1.03 


O TTT/ \ 

(J3 lll(l*) 


1 0£? 

1.36 


-5.8 


1 O TA 

18.79 


C A 

-6.0 


4.667 


1 A O 

-10.2 


1 A1 
1U1 


1 /i ao 
14.02 


1 AO 

1.08 


U5.5 V((t)) 


1 A A 

1.40 


-5.1 


1 A 1 O 

19.12 


-5.1 


/I c £ o 

4.568 


7 
-8.7 


36 


14.52 


0.99 


06 V 


1.31 


-5.1 


19.62 


-4.2 


4.583 


-7.9 


40 


13.33 


1.05 


03 V 


1.37 


-5.4 


18.73 


-5.7 


4.667 


-9.9 


Sn53 


14.47 


1 10 
1.12 


Oo.5 V 


1 A O 

1.43 


A A 

-4.9 


1 A 07 

19.37 


/i O 

-4.8 


a e 1 E 

4.515 


O A 

-8.0 


104 


13.02 


1.06 


03 111(1) 


1.39 


-5.8 


18.82 


-6.1 


4.667 


-10.3 


A o 

A2 


12.53 


1.04 


03 V 


1.36 


-5.4 


17.93 


-6.4 


4.667 


-10.7 


A3 


1 aa 
13.09 


1 A A 
1.04 


U3 lll(l^) 


1 07 

1.37 


K O 

-5.0 


1 O OA 

18.89 


K A 

-5.9 


4.667 


1 A 1 

-10.1 


I IT 

II ( 


1 /i 1 7 
14.17 


1 1 A 
1.10 


U6 V 


1 A O 

1.42 


-5.1 


1 A 07 

19.2/ 


-5.1 


/I coo 

4.583 


O 7 

-8.7 


Sn25 


12.23 


1 A O 

1.42 


bl lab 


1 £ A 

1.60 


-0.0 


1 O 70 

18. 16 


7 O 

-7.8 


A O /) O 

4.342 


A 7 

-9.7 


Sh64 


13.58 


1.15 


03 V((f)) 


1.47 


-5.4 


18.98 


-5.9 


4.667 


-10.1 


16 


13.53 


1.08 


03 V 


1.40 


-5.4 


18.93 


-5.6 


4.667 


-9.8 


102 


15.32 


1.10 


08.5 V 


1.41 


-4.9 


20.22 


-3.9 


4.515 


-7.0 


Sh57 


13.28 


1.06 


03 111(f) 


1.39 


-5.8 


19.08 


-5.8 


4.667 


-10.0 


108 


13.71 


1.10 


05.5 V 


1.42 


-5.2 


18.91 


-5.5 


4.597 


-9.2 


Shl8 


12.65 


1.21 


03.5 If 


1.53 


-6.3 


18.95 


-7.1 


4.597 


-10.8 


Sh24 


14.27 


1.11 


06 V 


1.43 


-5.1 


19.37 


-5.0 


4.583 


-8.6 


Sh49 


14.67 


1.10 


07.5 V 


1.42 


-5.0 


19.67 


-4.6 


4.541 


-7.9 


Sh47 


12.72 


1.14 


04 V 


1.47 


-5.8 


18.52 


-6.7 


4.643 


-10.8 


Sh23 


12.70 


1.14 


OC9.7 la 


1.39 


-6.0 


18.70 


-6.4 


4.481 


-9.3 


Sh22 


13.21 


1.06 


03 111(f) 


1.39 


-5.8 


19.01 


-5.9 


4.667 


-10.1 


Sh21 


14.75 


1.12 


06 V((f)) 


1.44 


-5.1 


19.85 


-4.6 


4.583 


-8.2 


Shl9 


13.65 


1.02 


03 V((f)) 


1.35 


-5.4 


19.05 


-5.3 


4.667 


-9.5 



a From spectral type and luminosity class. 

b Computed using our apparent spectroscopic distance modulus corrected by 4.3 x [E(B — V) — 1.39]. 
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Table 4. Stellar Distances to NGC 3603 



Study 


Method 


E(B - V) 


Adopted Ry 


Apparent DM a 


Distance 














Shcr 1965 


Main-sequence fitting 


1.42 


3.0 


17.0 


3.5 


Moffat 1974 


Main-sequence fitting 


1.32 


3.1 


18.7 


8.1 


Melnick & Grosbol 1982 


Main-sequence fitting 


1.38 


3.1 


17.9 


5.3 


Moffat 1983 


Spectroscopic parallax 


1.44 


3.2 


18.8 


7.0 ±0.5 


Melnick et al 1989 


Main-sequence fitting 


1.44 


3.2 


18.9 


7.2 


Crowther & Dcssart 1998 


Spectroscopic parallax 


1.23 


3.2 


18.8 


10.1 


Pandey et al. 2000 


Main-sequence fitting 


1.48 


3.1/4.3 b 


19.0 


6.3 ±0.6 


Sagar et al. 2001 


Spectroscopic parallax 


1.44 


3.1 


18.8 


7.2 ±1.2 


Sung & Bessell 2004 


Main-sequence fitting 


1.25 c 


3.55 


18.6 


6.9 ±0.6 


This Study 


Spectroscopic parallax 


1.39 


3.1/4.3 b 


19.1 


7.6 



a Apparent distance modulus computed using the quoted true distance modulus and the reddening correction 
made in each study. 

b A normal Ry of 3.1 is applied for the foreground reddening E(B — V) = 1.1, with a value of Ry applied to 
extinction within the cluster, i.e., A v = 3.1 x 1.1 + [(E(B - V) - 1.1) x 4.3] 

Tor the central cluster. 



Table 5. Kinematics Distances to NGC 3603 



Study Distance (kpc) 

Goss & Radhakrishnan 1969 7.1 a 

van den Bcrgh 1978 6.8 ± 0.9 a 

De Pree et al. 1999 7.0 b 

Russeil 2003 7.9 



"■Corrected to Ro = 8.5 kpc. 

b De Pree et al. 1999 quote a value of 6.1, but 
one of the co-authors, Miller Goss (private com- 
munication), writes that the value was a mistake, 
and recomputes a value of 7.0 kpc. 



